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Abstract. 

We study the interaction of an equal mass binary with an isothermal circumbinary 
disk motivated by the theoretical and observational evidence of the formation of 
massive black holes binaries surrounded by gas, after a major merger of gas-rich 
galaxies. We focus on the torques that the binary produces on the disk and how 
the exchange of angular momentum can drive the formation of a gap on it. We 
propose that the angular momentum exchange between the binary and the disk is 
through the gravitational interaction of the binary and a (tidally formed) global non- 
axisymmetric perturbation in the disk. Using this interaction, we derive an analytic 
criterion for the formation of a gap in the disk that can be expressed either via 
the characteristic velocities of the binary-disk system or in terms of the structural 
parameters structural parameters h/a and M(< r)/Mbi n . Using SPH simulations we 
show that the simulations where the binary opens a gap in the disk and the simulations 
where the disk does not have a gap are distributed in two well separate regions. Our 
analytic gap-opening criterion predicts a shape of the threshold between this two 
regions that is consistent with our simulations and the other ones in the literature. 
We propose an analogy between the regime without (with) a gap in the disk and the 
Type I (Type II) migration that is observed in simulations of planet-disk interaction 
(binaries with extreme mass ratios), emphasizing that the interaction that drives the 
formation of a gap on the disk is different in the regime that we analyze (comparable 
mass binary). 
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1. Introduction 

It's well known that galaxies are systems that can strongly interact gravitationally with 
each other and even in some cases merge between them. Also, it's widely accepted the 
fact that galaxies with a significant bulge host a massive black hole (MBH) at its center 
(Richstone et al 1998). Therefore, is natural to conclude that these MBHs will interact 
after a major merger (Milosavljevic & Merritt 2001, Milosavljevic & Merritt 2003) and 
possibly collide. The possible coalescence of MBHs was first considered by Begelman et 
al (1980) in a study of the long-term evolution of a black hole binary at the center of a 
dense stellar system. Initially, dynamical friction brings the two black holes towards the 
center of the system and the resulting MBH binary continues to shrink via three-body 
interactions with the surrounding stars. This three-body interactions tends to eject stars 
from the central region causing the merger eventually to stall ("last parsec" problem) 
unless some additional mechanism is able to either extract angular momentum from the 
MBH binary or refill the stellar "loss-cone" (Khan et al 2011). 

There are theoretical (Barnes & Hernquist 1992, 1996; Mihos & Hernquist 1996; 
Barnes 2002; Mayer et al 2007, 2010) and observational (Sanders & Mirabel 1996; 
Downes & Solomon 1998) evidence that in the mergers of gas-rich galaxies a large 
amount of gas can reach the central regions of the newly formed system. This gas 
induces the formation of a MBH binary (Kazantzidis et al 2005, Mayer et al 2007, 
Chapon et al 2011) and can drive the final coalescence of the binary (Escala et al 2005; 
Dotti et al 2006; Cuadra et al 2009; Rodig et al 2011; Escala & del Valle 2012). 

In all these studies that follow the evolution of a binary embedded in a gas 
environment is assumed that the gas lies in a disk (circumbinary disk) and the binary 
has a non-extreme mass ratio (q ~ 1). However, they find very different times scales 
of coalescence for the binary because they explore very different disk to binary mass 
ratios (Mdi S k/ Mu n ). When the mass of the disk is much greater than the binary mass 
it's found that the coalescence time of the binary is on the order of a few initial orbital 
times (Escala et al 2005; Dotti et al 2006). On the other hand, for disk's with negligible 
masses, compared to the binary mass, the coalescence time is on the order of several 
thousand of local orbital times (Artymowicz & Lubow 1994; Ivanov, Papaloizou, & 
Polnarev 1999; Armitage & Natarajan 2002; Milosavljevic & Phinney 2005), which for 
Mbh > 1O 7 M is even longer than the Hubble time (Cuadra et al 2009). 

These two regimes of fast (few orbital times) and slow (several thousands of orbital 
times) migration are also observed in simulations of planets embedded in circumstellar 
disks (eg. extreme mass ratios; Ward et al 1989, Ward 1997, Bate et al 2003, Armitage 
& Rice 2005; Baruteau & Masset 2012, Kocsis et al 2012a, 2012b). In simulations of 
planet migration in protoplanetary disks (q « 1), these two regimes are defined as 
Type I and Type II migration; In the Type I regime the perturbation induced by the 
planet in the circumstellar disk remains small and permits a fast migration of the planet 
(of the order of a few orbital times) with a characteristic time scale that is inversely 
proportional to the planet's mass {t m ig ration oc M" 1 ). When the Hill radius of the planet 
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is greater than the local pressure scale height h (Rmn » h) the perturbation induced 
by the planet in the disk become important and a gap begins to form, leading to coupled 
evolution of the planet and the disk on a viscous time scale and therefore, a much longer 
migration time (Type II migration). In this paper we'll extend the same terminology to 
the case of non-extreme mass ratios binaries (q ~ 1) interacting with a disk by referring 
to the fast (slow) coalescence of binary as Type I (II) migration. 

For a comparable mass binary in a disk, as in the star-planet-disk system, the 
threshold between the Type I and Type II migrations is determined by the formation of 
a gap in the disk. Our aim in this paper is establish a gap opening criterion that allows 
us to differentiate these two migrations regimes and determine in which binary-disk 
systems is expected an efficient coalescence of the MBH binary. For that purpose, we 
will use a simple model for the binary-disk interaction to derive an analytic expression 
for the gap-opening criterion. We will test this criterion against a variety of simulations 
with different model parameters for the binary to disk mass ratio and thickness of the 
disk. 

This paper is organized as follows. After the derivation of an analytic gap-opening 
criterion in §2, we present in §3 the numerical method and set of numerical simulations 
that we use to test our analytic criterion. In §4 we describe how we identified the 
formation of a gap in our simulations and in §4.1, we test our analytic gap opening 
criteria against the numerical simulations. All the results that we find are discussed in 
§5 and in this section are also compare our results with the results of other authors in 
the literature. Also, in §5 we discuss the implications of our results in real astrophysical 
systems. Finally, our conclusions are presented in §6. 

2. ANALYTIC GAP-OPENING CRITERIA 

The interaction of a binary with a circumbinary disk in cases where the total gas mass 
is typically much smaller than the mass of the primary, has been widely studied in 
the context of star/planet formation (Lin & Papaloizou 1979; Goldreich & Tremaine 
1982; Takeuchi & Lin, 2002; Armitage & Rice 2005; Baruteau & Masset 2012). In these 
studies is used a linear approximation for the equations of motion of the system and it's 
found that the interaction between the secondary and the disk is controlled by the sum 
of the torques arising from the inner and outer Lindblad and corotation resonances. 
This approach applies to binaries where the primary is much more massive than the 
secondary (q = M 2 /M\ <C 1) and leads to predictions of the gap structure that are 
consistent with simulations within the same regime (Ivanov et al . 1999;Bryden et al 
1999 ;Armitage & Natarajan 2002; Nelson & Papaloizou 2003; Haiman et al . 2009; 
Baruteau & Masset 2012). Because of the success of this analysis in the planetary regime 
(q <C 1), it has been extrapolated to other cases where q ~ 1 (Artymowicz & Lubow 
1994, 1996; Gunther & Kley 2002; MacFadyne & Milosavljevic 2008) without consider 
the global non-linear perturbation that is produced by the binary gravitational field in 
the regime q ~ 1 (Shi et al 2012). This non-linear perturbation breaks the validity of 
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the linearization of the equation of motion. 

In this paper we study the binary-disk interaction in the regime q ~ 1 without any 
assumption of linearity. For this purpose we explore if it's possible, due to the tidal 
nature of the binary-disk interaction, that the gap-opening process can be described 
by the interaction of the binary with a strong non-axisymmetric perturbation on the 
disk instead of a resonant process that appears in the linear approximation. This 
type of interaction between a strong non-axisymmetric perturbation and a binary was 
investigated in Escala et al (2004, 2005). They show that, when the gravitational 
influence spheres (i.e. Hill spheres) of the two component of the binary overlap, a 
strong tidal non-axisymmetric perturbation is produced in the disk with an ellipsoidal 
geometry for an equal mass binary (q — 1). They also found that the symmetry axis 
of this strong non-axisymmetric perturbation is not coincident with the binary axis but 
lags behind it, producing a gravitational torque on the binary that is responsible for an 
angular momentum transport from the binary to the disk. 

In this paper we will restrict to the case of an equal mass binary (g = 1) interacting 
with a disk, leaving the case q ~ 1 but q ^ 1 for a companion paper. In the case q — 1, is 
justified to assume that the strong non-axisymmetric perturbation produced in the disk 
has an ellipsoidal geometry. Therefore, we study the gravitational interaction between 
the binary and an ellipsoid to derive analytically the exchange of angular momentum 
between the disk and the binary. This binary-ellipsoid system can be treated as an 
equivalent one body problem subject to an external gravitational potential caused by 
the ellipsoidal perturbation. The gravitational potential of a uniform ellipsoid is given 
by 

y, z) = it G p ( a x 2 + /3 y 2 + ^z 2 + Xo) , (!) 

with ao, Po, 70 and xo constant given in Lamb (1879) which depend on the ratio between 
the principal axes of the ellipsoid. For the symmetry of the ellipsoidal perturbation 
observed in Escala et al (2004, 2005) 70 = /3q and xo — 0. Therefore, the Lagrangian of 
the binary-ellipsoid system in cylindrical coordinates has the form 

2 / , , \ 2 / , \ 2" 
+ 



L = -a 
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2 Gfip(a r 2 cos 2 (A0) + f3 r 2 sin 2 (A0) + (3 Q z 2 ) , (2) 

with pL = MiM 2 /M b in being the reduced mass of the binary and A0 the angle between 
the major axis of the ellipsoid and the binary axis. This angle is approximately constant 
over the evolution of the system (Escala et al 2004, 2005). From this Lagrangian we 
derive the Euler-Lagrange equation for the coordinate 0, 



dt dt \ dt 



_ nr 2 ) P G(f3 - a ) cos(A0) sin(A0) , (3) 
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where we use \i = M bin /4. This equation expresses the decrease of the angular 
momentum of the binary /^n and therefore, the injection of angular momentum (torque) 
from the binary to the disk. 

To derive a criterion for the opening of a gap in the disk we compare the gap- 
opening time scale (determined by the torque that the binary exchange over the disk) 
with the gap-closing time scales (Goldreich & Tremaine 1980). The angular momentum 
AL that must be added to the gas to open a gap of radius Ar in a disk with thickness 
h, is of the order of AL p (Ar) 2 hrv. The torque that the binary exchange over the 
disk is r = — dl bin /dt = r 2 p (M bin /4) G ir (j3 — a ) cos(A0) sin(A0). This torque injects 
an angular momentum AL on a time scale At opcil = AL/t and therefore, this is the 
characteristic time scale to open a gap. This tendency is opposed by viscous diffusion, 
which fills up a gap of with Ar on a timescale At c i osc = (Ar) 2 jv (Goldreich & Tremaine 
1980) where v is the turbulent viscosity of the gas that can be parametrised assuming 
the standard a-prescription v = a ss c s h of Shakura & Sunyaev (1973) where a ss is the 
dimensionless viscosity parameter and c s is the sound speed of the gas. Gap formation 
occurs when A£ opcn < At close . We also assume that c s /v ~ h/r with v the circular 
velocity of the binary-ellipsoid system. Here, based on previous numerical results (Escala 
et al 2004, 2005), we assume that the binary and the ellipsoid (strong non-axisymmetric 
perturbation) co-rotate in a circular orbit. If we define t> bin = (GM hin / (Sr)) 1 ^ 2 is 
straightforward to find that the binary will open a gap in the disk if 



where / = f(A(f), a , (3 , a ss ) = 27r(/5 — cxo) cos(A0) sin(A0)o; ss 1 is a dimensionless 
function of the geometric parameters of the ellipsoid (etc A)), the offset angle A<p and the 
viscosity parameter a ss . The geometric parameters of the ellipsoid are approximately 
constant due to the self-similar behavior of the ellipsoid (Escala et al 2004, 2005). This 
formulation of the gap opening criterion contains information of the relative weight of 
the characteristic rotational speed of the binary-disk system (v), with the characteristic 
rotational speed of the isolated binary (r> b m) and the thermal state of the disk (c s ). 

Regardless of the exact geometry of the strong non-axisymmetric density 
perturbation, the gravitational torque from the perturbed background medium onto 
the binary will have the form r = r 2 pGfiK, where K is a parameter that depends 
on the geometry of the density perturbation. Therefore, the dimensionless function / 
has in general the form / = 2K/a ss that represents the relative strength between the 
gravitational and viscous torque. 

If we assume that the binary-disk system is in rotational support against the overall 
gravitational potential on a circular orbit, we can also express this same criterion in 
terms of the disk structural parameters and the mass of the binary-disk system: 




(4) 




(5) 



where M(< r) is the total mass enclosed by the orbit of the binary. 
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3. INITIAL CONDITIONS AND NUMERICAL METHOD 

In this section we present simulations that follow the evolution of a binary embedded in a 
gaseous disk. We use a natural units system in our simulations: [mass] = l, [distance] = 1 
and we set the gravitational constant G=l. In our internal units, the initial radius of 
the disk is -Rdisk = 30, and the total gas mass is M disk = 33. The density of the disk is 
given by 



where R c = 3, and h c are the radius and thickness of the central zone of the disk 
where the surface density is constant. With this density we obtain a surface density 
E(r) = constant for R < R c and S(r) oc R" 1 if R > R c . The vertical distribution of 
the disk changes over the evolution of the system, but the initial vertical distribution 
chosen at least prevents an initial vertical collapse on the disk. The binary system have 
an initial circular orbit of radius r , coplanar and co- rotating with the disk, a total mass 
M bin , and a mass ratio between the two components of the binary q = M 2 /M 1 = 1. 
We explore the parameter space r G [0.5,4], M h - m G [1,33], and h c G [0.8,3] with 25 
simulations (see table 1). 

In addition to the disk and the binary we include a fixed Plummer potential 
(Plummer 1911) that helps to stabilize the disk and also, when we apply our result 
on the study of SMBH binaries in section §5, will mimic the existence of an external 
stellar component. The total mass of the Plummer potential is M p = 6.6 ( ~ %12 of 
the total mass of the disk) and its scale radius is R p = 19.13. The applicability of the 
gap-opening criterion derived in §2, remains valid even in the presence of this Plummer 
potential because its mass can be included in the total mass enclosed by the orbit of 
the binary M(< r). 

We model the gaseous disk with a collection of 2 x 10 5 SPH particles of gravitational 
softening 0.1. We use a stable (Q > 1) isothermal disk to avoid fragmentation, to 
simplify the testing of our analytic criterion derived in the section §2. In table 1 we 
tabulate the minimum value of the parameter Q of Toomre of each simulation. Since the 
total potential of the system (Binary-Disk-Plummer) is non-axisymmetric and therefore 
lacks of a well defined circular velocity, we assume initially a symmetric representation 
of the binary potential to compute the circular velocity of the disk (see Appendix A for 
details). For the binary we use 2 collisionless particles with gravitational softening of 
0.1. We run each simulation for 10-15 binary orbital times using the SPH code called 
GADGET2 (Springel 2005, Springel et al 2001). 

In addition, to test the numerical convergence of the results, we run simulations 
with 1 million SPH particles and find that the conclusions found in §4 for our analytic 
gap-opening criterion do not change (Appendix B). 




(6) 



(7) 



Gap-Opening criteria 



7 



Table 1: Run Parameters 



RUN 


r o/Rdisk 


M(< r)/M bin 


(Cs/Vbin) 2 


Qmin 


Al 


0.133 


0.05 


0.088 


2.62 


Bl 


0.1 


0.06 


0.080 


2.50 


B2 


0.1 


0.06 


0.159 


2.50 


B3 


0.1 


0.06 


0.199 


2.49 


B4 


0.1 


0.06 


0.279 


2.48 


B5 


0.067 


0.06 


0.053 


2.50 


CI 


0.1 


0.1 


0.133 


2.22 


C2 


0.067 


0.1 


0.089 


2.20 


C3 


0.1 


0.1 


0.199 


2.20 


C4 


0.067 


0.1 


0.178 


2.20 


C5 


0.067 


0.1 


0.222 


2.20 


C6 


0.1 


0.1 


0.465 


2.19 


C7 


0.067 


0.1 


0.311 


2.19 


Dl 


0.1 


0.12 


0.159 


2.14 


D2 


0.05 


0.12 


0.080 


2.15 


D3 


0.05 


0.12 


0.199 


2.14 


El 


0.067 


0.2 


0.178 


2.01 


E2 


0.033 


0.2 


0.088 


2.01 


E3 


0.033 


0.2 


0.219 


1.99 


E4 


0.067 


0.2 


0.623 


1.98 


Fl 


0.017 


0.4 


0.090 


1.85 


F2 


0.033 


0.4 


0.175 


1.85 


F3 


0.067 


0.4 


0.356 


1.85 


Gl 


0.017 


2.0 


0.452 


1.70 


G2 


0.067 


2.0 


1.779 


1.70 



4. RESULTS 

In order to test our analytic gap opening criterion against SPH simulations, we first 
need a criterion that determines numerically which of our simulations opens a gap in 
the disk and which ones don't. To determine whether a gap is opened or not, we explore 
the radial density profile of each simulation and the evolution of the binary separation. 

If the binary opens a gap in the disk, it must have a flow of gas from the central to 
the outer regions of the disk. This flow produces a "pile-up" of gas on the perimeter of 
the gap. This type of pile-up is also expected to form on the evolution of a circumbinary 
disk that interacts with an unequal mass binary but in this regime (q <C 1) is associated 
to a drop in the density of the disk and not necessarily to the formation of a gap (Kocsis 
et al 2012). The pile-up formed in the evolution of our simulations (q = 1) is represented 
by a peak in the density profile and its maximum is correlated to the existence of a gap 
in the disk. In the figured] we show that the density plot has a peak or pile- up if the disk 
has a gap. For the cases that the disk doesn't have a gap, this pile-up is less prominent 
or doesn't exist. 
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Figure 1. In the top: the density projection in the plane x-y of the simulation C2 
(left) and F3 (right). In the bottom: The density profiles of the simulation C2 (left) 
and F3 (right). In this set of graphics is shown the correlation between the simulations 
with a gap (without a gap) and a density profile with a pile up of gas whose maximum 
value is p pea k > 0.015 (/9 pea k < 0.015) (see section §4) 



We also explore if the gap can be identified as a drop in the azimuthal mean density. 
Although in some cases the existence of a gap is correlated with a drop in the mean 
azimuthal density, in many other cases there are streams of gas or inner disks around 
individual components of the binary that hide the existence of a gap in an azimuthally 
averaged density plot. Therefore, we decided to use the peak in the density profile as 
our first indicator for the existence of a gap. In addition to this, we follow the evolution 
of the binary separation. This separations tell us whether the binary is migrating or not 
to the inner regions of the disk on a time scale comparable to its orbital time t or b- In 
order to the binary migrates is necessary an efficient angular momentum transfer that 
cannot be reached if the material that participates in this angular momentum exchange 
is diffused to the outer regions as in the open-gap case. Therefore, a binary that is 
migrating to the center of the disk on a time scale ~ t or b, is only consistent with a disk 
that does not have a gap (Type I migration). 

Considering the density profile and the possible migration of the binary to the inner 
region of the disk, we called opened to a simulation in a given time t, if it has a peak or 
pill-up in density with p pea k > 0.015 (in internal units) and if the binary separation do 
not decreases by more than 10 % in a orbit. If a simulation in a time t have a maximum 
density of the pile up p pea k < 0.015 and the binary separation decreases more than 10 
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% in a orbit, we called closed. Otherwise we called it indefinite. 

We analyze our simulations at the times t in which the binary complete 2, 3, 5, 7, 
10 and 15 orbits. For all the simulations in each of these time intervals, we define it as 
opened, closed or indefinite if the gap in the disk is open, close or indefinite respectively. 

4-1. Testing the Gap-opening Criterion 

In order to test our analytic gap-opening criterion we plot the velocity parameters of 
the system ( (t>bi n / v) 2 , ( c s / v) 3 ). This parameters are present in one of the two 
equivalents gap-opening criteria (equation HJ). In a first attempt we try to determine 
the values of the parameters ( 1 + M(< r) / M^- in )~ 1 and (h / r) 3 that are present in the 
other gap-opening criterion (equation [5]) but were poorly defined because the structural 
parameters of the disk (thick, radius and density profile) can vary strongly (variations 
on the order of 30-70% ) in the vicinity of the binary, in time-scales shorter than the 
orbital time-scale of the binary, due to the local action of the gravitational potential of 
the components of the binary. 

The velocity parameters are derived from the dynamical state of the binary that, 
unlike the structural parameters of the disk, have variations in times scales greater 
or at least comparable, to the time scale of the orbit of the binary. Therefore, to 
compute the structural parameters of the system we derive them from the velocity 
parameters. These structural and velocity parameters are related through the relations 
8(M(< r)/M bin ) + 1 = {v/v h - m ) 2 and {h/rf = (c s /vf. 

In figure [21 we plot two groups of simulations; the opened and the closed simulations 
according to our definition. The opened simulations are represented by red open circles 
and the closed simulations are represented by blue filled circles. We find that the 
opened and closed simulations are distributed in two different regions. The interface 
that separates the set of simulations that are opened from the set of simulations that 
are closed can be modeled by a straight line (which slope is f as we will show in the next 
section). The linear nature of the interface indicates that the functional dependence of 
our gap-opening criterion on the parameters of the disk-binary system works fine, in the 
range of parameters that we explore. This encourages us to explore in more detail the 
significance of this linear interface. 

An important consequent of Figure 2 is that the extension of the standard gap- 
opening criterion derived for the planetary regime (q <C 1) to the regime q ~ 1 doesn't 
predicts the right shape of the interface line between the opened and closed simulations. 
For example the standard gap opening criterion in Lin & Papaloizou (1986) (Equation 
(4)) for the case q = 1 can be expressed as (h/r) 3 < (l/40a ss ) 3//5 which corresponds to 
an horizontal line in Figure 2, which clearly cannot explain the distribution of opened 
and closed simulations in [2j 




Figure 2. The figure shows the cubic ratio between the sound speed of the gas and 
the rotational velocity of the binary-disk system (c s /v) 3 plotted against the quadratic 
ratio between the rotational velocity of the isolated binary and the rotational velocity 
of the binary-disk system (v hin /v) 2 . As (vbin/v) 2 = (1 + 8M(< r)/Mbin) -1 and 
(c s /v) 3 = (h/r) 3 we also include this axis labels where M(< r)/Afbm is the enclosed to 
binary mass ration and h/r is the ratio between the disk thickness and half the binary 
separation (r — a/2). The red circles are simulations where the binary has opened a 
gap in the disk (opened simulations) and the blue filled circles are simulations where 
the disk does not have a gap (closed simulations) (see section §4). The green line is the 
threshold between the opened simulations and the closed simulations that is predicted 
by our analytic gap-opening criterion. Below the green line are the opened simulations 
and above the line are the closed simulations. The slope of the interface is the function 
/(A</>, a o ,0o,a ss ). 



4-2. Determining an average value for f(A<p, cto, /3q, «ss) 

The interface-line that was shown in the previous section separates the parameter space, 
between the opened and closed simulations and can be interpreted as the critical case 
in which a simulation have equals opening and closing times (i.e At open = At c i ose ). The 
slope of the interface-line that separates the set of opened simulations from the set of 
closed simulations in the velocity variables graph (figure [2]) is 



c s \ 3 ( h 



3 



M(< r) 

1 + ° 



bin 



(8) 



We find a value for the slope of the interface-line approximately of m = 0.095. 

From the equation 0] for the case At open = At c i ose we can estimate that the 
dimensionless function / that contains the information of the relative strength between 
the gravitational and viscous torques , is on average /(A(/>, a , (3 , a ss ) = m. Replacing 
this numerical value for / in the equations H] and [5] we can express the gap-opening 
criterion for an equal mass binary as 



Gap-Opening criteria 



11 



Vbi: 



- < 0.095. 



or 



7;' 



3 r 



M(< r) 



M, 



bin 



< 0.095 



(9) 



(10) 



As / oc a" 1 any changes on the viscosity parameter a ss will change the slope of 
the interface-line of the figure [2j If we increase the value of a ss the slope will be less 
steep and the number of closed simulations will increase. This is consistent with the 
fact that with an increase of the viscosity it will be harder for the binary to open a gap 
on disk. In our simulations we can estimate a ss from the value of the SPH parameter of 
artificial viscosity a spn (Artymowicz & Lubow 1994, Murray 1996, Lodato & Price 2010, 
Taylor Sz Miller 2011). The value that we estimate ranges between a ss ~ 0.008 ~ 0.016. 
However, independent to the exact value; the functional dependence of our gap-opening 
criterion remains unchanged. 



4-3. Transition from closed regime to opened regime 

Only one of our 25 simulations (simulation C5 in table [3]) evolves from a state where 
the disk has no gap (closed) to a state where the disk has a gap (opened). The binary 
in this simulation initially migrates to the center of the disk driven by the action of 
the tidal torque. As the binary migrates to the center, eventually reaches a separation 
where the torque exerted into the disk is enough to open a gap. It is important to 
emphasize that this is only a particular case that already start near the opened region 
and doesn't mean that any binary as migrates to the center will be able to open a 
gap. Moreover, in most cases studied, we have binaries that as migrate to the center 
never fulfil the condition to open a gap, because the mass ratio Mu n /M(< r) doesn't 
grows fast enough compared to the growth of the thickness h/r. It's interesting to 
visualize the transition of the simulation C5 between the opened and closed gap regions 
in the velocity parameters space ((vun/v) 2 , (c s /v) 3 ) (figure [3]). In this plot we find that 
the three regimes of evolution in this simulation are distributed consistently with our 
criterion derived in §2, with the closed points above the threshold line (/ = 0.095), the 
opened points below the threshold and the indefinite ones, oscillating around it. 




Figure 3. The figure shows the cubic ratio between the sound speed of the gas and the 
rotational velocity of the binary-disk system (c s /v) 3 plotted against the quadratic ratio 
between the rotational velocity of the isolated binary and the rotational velocity of the 
binary-disk system (vbin/^) 2 - We plot three groups of points representing three groups 
of different regimes for a same simulation in different times (C5). The three regimes are 
opened (red open circles), closed (blue filled circles) and a transitory indefinite regime 
(green triangles) according to our definition (see section §4). The green line is the 
interface between the opened population of parameters (below the line) and the closed 
population of parameters (above the line). The slope of the interface is the function 
/(A</>, ceo,0Q, a ss ). The three regimes are distributed in the plot consistently with the 
criterion. 



5. DISCUSSION 

In our set of simulations we identify two regimes of binary-disk interactions: with a gap 
(opened) and without a gap (closed). We found that in general less massive and thinner 
disks are strongly perturbed by the binary and a gap is formed in them. On the other 
hand, an increase of the mass or thickness of the disk tends to close this gap and in 
some cases, even preclude its formation. We also found that the possibility of formation 
of such a gap is more sensitive to the thickness than to the mass of the disk. This is 
in agreement with our analytic criterion that has a linear dependence in the mass ratio 
M bin /M(< r) and a cubic dependence in the ratio between the pressure scale height of 
the disk and the binary separation (h/a). 

We find good agreement between our test-simulations and our analytic criterion, 
this support the role of the interaction between the binary and a strong non- 
axisymmetric perturbation as the responsible of the formation of the gap in the disk, 
instead of a resonant process. However, the resonances that operates in regions far from 
the binary where the linearization of it gravitational potential is a good approximation 
(r >> a), play an important role in the later evolution of the gap, "clearing" it edges 
and driven the angular momentum transport from the binary to the disk when the gap is 
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already formed (Artymowicz & Lubow 1994; MacFadyen & Milosavljevic 2008; Cuadra 
et al 2009; Roedig et al 2012). 

To improve the robustness of our study in figure H] we plot the parameters 
M(< r) I Mbi n and h / r (being r = a/2 with a the binary separation) for our simulations 
and the close and opened simulations from the literature (Artymowicz & Lubow 1994; 
Artymowicz Sz Lubow 1996; Escala et al 2005; MacFadyen & Milosavljevic 2008; Cuadra 
et al 2009). In this figure we also plot the interface-curve (At opcn = At closc ) that for this 
space of parameters has an inverse-cubic form. For the cases of simulations with non-self 
gravitating disks (Artymowicz & Lubow 1994; Artymowicz & Lubow 1996; MacFadyen 
& Milosavljevic 2008) we assume M(< r)/Mu n = 0. All the opened simulations that 
we find in the literature begins with a gap in the disk (M(< ro)/Mbi n = 0) or the disk 
have a low mass compare to the binary mass, but we dont find any simulation where 
the formation of a gap is study for a massive disks (M(< ro)/Mu n > 1). In the figure 
[4] we can see that the simulations in which the disk begins with a gap, or it's open in 
the evolution of the simulation, (red open points I, II and III) are below the interface 
curve, consistent with the gap-opening criterion. The closed simulations (squares in 
figure Hj) represent the final evolution of simulation that begins with initial condition 
as the used by Escala et al 2005 and Dotti et al 2006. These simulations are above 
the interface line and are also consistent with the gap opening criterion. Therefore all 
the simulations from papers that we refer from the literature are consistent with the 
gap-opening criterion independently of the different equations of state that are used for 
the disks. 

In the final evolution of the simulations from Escala et al (2005), squares in figure 
HI is found that although the binary shrinks their separation up to distances where its 
gravitational potential begins to dominates locally and the enclosed mass is comparable 
or smaller than the binary mass (M(< r)/Mb; n < 1), the binary fails to form a gap in 
the disk. This is because the scale height of the disk remains roughly constant as their 
separation a (equals to 2r) decreases, resulting in an increases of h/a that compensates 
the decrease in M(< r)/Mbi n in such a way that the system always lies above the 
threshold line / = 0.095. 

We expect that a situation like the one found in Escala et al (2005), should happen 
in galaxies with a gas rich nucleus, such as nuclear disks in ULIRGs (Downes & Solomon 
1998), Submillimiter Galaxies (Chapman et al 2003; Chapman et al 2005; Tacconi et al 
2006; Swinbank et al 2010) and in general, on protogalaxies at the early universe. This is 
because both observations (Downes & Solomon 1998; Genzel et al 1998) and simulations 
(Mayer et al 2010; Bournaud et al 2011) found typical disk's scale height larger than 
several tenths of parsecs and gaseous masses larger than 10 9 M Q (more than the mass 
of most massive MBHs binaries). Moreover, we determine f in our gap opening criteria 
(equation HI) with simulations using an a ss ~ 0.01 (f oc a" 1 ) and gas in massive nuclear 
disks is expected to be globally unstable. In such a case, the torques are significantly 
larger, with a ss of order unity (Krumholz et al 2007; Escala 2007). Therefore, in this 
situation we expect that to open a gap will be even harder than what we found in our 
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• Close 
Open 
— f= 0.095 




Figure 4. The ratio between disk thickness and binary separation (h/a) plotted 
against the enclosed to binary mass ratio (M(< r)/Mbi n ), for the simulations with a 
gap in the disk (open red circles) and without gap in the disk (blue filled circles). The 
continuous green curve represent the interface of this two populations of parameters 
(f = 0.095). The other points in the figure, are simulations already in the literature: 
Escala et al 2005 (1), Artymowicz & Lubow 1994 (I), Artymowicz & Lubow 1994 and 
MacFadyen & Milosavljevic 2008 (II), Cuadra et al 2009 (III). The points with arabic 
number labels, are simulations where the disk don't have a gap and the points with 
roman number labels are simulations where the disk have a gap. All the points are 
consistent with our analytic gap-opening criterion. The qualitative behavior predicted 
by the opening criterion is reproduced in this graph, as h/a decreases, the binary-disk 
interaction tends to open a gap. Also if we decreases M(< r)/AIti n , the interaction 
tends to open a gap again. 

simulations. 

On the other hand, in mergers of gas poor galaxies we expect MBHs binaries with 
less massive and thinner circum-binary disks in which a gap will be easily opened. In 
such cases, the MBH-disk interaction will be better described by Type II simulations 
(Artymowicz & Lubow 1994; Artymowicz & Lubow 1996; MacFadyen & Milosavljevic 
2008; Cuadra et al 2009). However, in order to make concrete predictions, is needed 
high-resolution galaxy mergers simulations that includes feedback processes from star 
formation and BH accretion, to study under which conditions the criteria given by 
equation H] is fulfilled or not and therefore, the binary-disk interaction will be Type I or 
II. 

Our study has also implications for the formation of gaps around binary protostars, 
in which both Type I and II interactions can be present at different epochs of their 
formation. In simulations of early epochs in the formation of binary stars, for example 
Boss (1982), Bate & Bonnell (1997) that studies the evolution of binary proto-stellar 
seeds within a collapsing cloud, is found that the interaction of the binary with the 
ambient gas is through a non-axisymmetric overdensity like the one studied in this 
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paper as Type I. On the other hand, in studies of the final stages in the formation of 
binary stars, simulations of proto-planetary disks around binary stars recurrently found 
gaps that produce disk truncation at r t ~ 3/2a (Artymowicz & Lubow 1994; Artymowicz 
& Lubow 1996; Gunther & Kley 2002). Between both limiting epochs in the formation 
of binary stars, it must be a transition epoch that is given by the threshold condition 
studied in this paper (equation H]). When this transition occurs in binary formation, 
there must be implications in the final separations and masses of the binaries, since it 
corresponds to a dramatic transition in both the migration and accretion timescales. 

6. CONCLUSIONS 

We study the interaction of an equal mass binaries with an isothermal gas disk. We 
focus on the torques that the binary produces on the disk and how this exchange of 
angular momentum can drive the formation of a gap on it. 

We propose that the angular momentum exchange between the binary and the 
disk is through the gravitational interaction of the binary and a tidally formed strong 
non-axisymmetric perturbation in the disk. From this gravitational torque we derive an 
analytic criterion to determine if the binary will open a gap or not in the disk. 

The analytic gap opening criterion that we derive depends on two dimensionless 
parameters (fbinAO 2 and (c s /t>) 3 where v is the rotational velocity of the binary-disk 
system, vu n is the rotational velocity of the isolated binary and c s is the sound speed 
of the disk. We also show that, through a transformation of variables, the criterion 
can be rewritten into other two parameters; the first is a function of the enclosed to 
binary mass ratio (1 + M(< r)/M bin ) and the second is the cubic ratio between the disk 
thickness and the binary separation (h/r) 3 . 

Using SPH simulations we show that the simulation where the binary opened a gap 
in the disk (opened simulations), and the simulations where the disk does not have a gap 
(closed simulations), are distributed in two well separate regions in the parameter spaces 
( {vbin/ v ) 2 i ( c s/ v) 3 ). Our analytic gap-opening criterion predicts a linear shape of 
the threshold between this two regions and our SPH simulations allows us to determine 
the value of it slope. The value that we find for the slope is roughly 0.095. 

To increase even more the confidence in our results we test our analytic gap opening 
criterion against simulations of the literature and we find a good agreement between 
our gap opening criterion and this simulations. 

Our simulations assume only one type of initial conditions for the binary-disk system 
and we leave for another publication the study of the possible effects of different initial 
conditions on the results. However, we emphasize that for an isolated binary-disk system 
the gap opening/closing criterion depends on the local parameters at a given time (e.g. 
the velocity of the binary-ellipsoid system) but how this local parameters evolves through 
time depends on the initial conditions, for that reason, our plan in the future is to vary 
them. 

Finally, we discuss how our results can be applied to the study of the formation 
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of gaps on the circumbinary gas disks around binary MBH and protostars. We discuss 
the implications of the formation of a gap on the migration and accretion timescales of 
these systems. Despite the generality of the process that drives the formation of a gap, 
we remark that is mandatory to explore how more realistic simulations could affect the 
gap-opening criterion. 

Acknowledgements. L del V research was supported by CONICYT Chile (Grant 
D-21090518) and DAS Universidad de Chile. A.E. acknowledges partial support 
from the Center of Excellence in Astrophysics and Associated Technologies (PFB 06), 
FONDECYT Iniciacion Grant 11090216. The simulations were performed using the 
HPC clusters Markarian (FONDECYT 11090216) and Geryon (PFB 06). 



Gap-Opening criteria 



17 



APPENDIX A: Initial Conditions 

In our set of simulations, the binary produces a strong non-axisymmetric component to 
the total gravitational field. Therefore, for the SPH particles near the binary there is not 
a well defined circular velocity v(r). This is an inherent problem for the initialization 
of this type of simulations and we chose to solve this by approximating the binary 
potential by an spherical mass distribution with the same total mass, for computing the 
initial circular velocity. These spherical representations allow us to compute an initial 
circular velocity for all the particles of the disk, but the system will not start in a perfect 
equilibrium. 

We probe different spherical mass distributions to study the effect of the initial 
velocity on the evolution of the system and then we select the mass distribution that 
produce the minor overestimation or underestimation of the equilibrium velocity. The 
spherical mass distributions that we study are: (1) spherical shell (2) point mass 
distribution and (3) sphere of constant density. The spherical shell (1) produces a 
discontinuity on the velocity that generates a ring-shaped gap due to an overestimation 
of the velocity in the region where the gap forms. The point mass distribution (2) 
overestimates the velocity for all the particles in the region r < ro because in this 
region the real enclosed mass is less than the mass that is assume by the point mass 
distribution. This overestimation opens an artificial gap in the center of the disk. The 
sphere of constant density (3) also generates an overestimation of the velocity in the 
region r < r but this overestimation is less intense than the overestimation due to 
the point mass distribution, because, for the homogeneous sphere, the enclosed mass 
decreases with radius. 

We select the sphere of constant density as the spherical mass distribution for 
computing the initial velocity, because approximates better the potential of the binary 
and the system starts closer to a rotationally supported equilibrium. The system is of 
course not in a perfect equilibrium and our approximation for the initial circular velocity 
leads to an initial readjustment of the density profile in about an orbital time. For that 
reason, we only analyze the simulations after 2 initial orbits of the binary, that is, after 
the system has relaxed into an quasi-equilibrium configuration. 

We also include a static Plummer potential to increase the stability of the disk. This 
static potential decrease the overestimation of the rotational velocity in the inner region 
of the disk reducing the initial loses of gas in this region and the artificial formation 
of a gap. In the edges of the disk the cutoff of the disk produces a pressure gradient 
that drives the expansion of the disk edges. This flow of gas is also reduced due to the 
gravitational influence of this static potential. We use an static Plummer potential that 
have a total mass of 12% the mass of the disk and its core radius is roughly at 0.6 -Rdisk- 
With this selection of parameters the Plummer potential's mass initially enclosed by 
the orbit of the binary is %50 the total enclosed mass. 
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APPENDIX B: Resolution Study 

Our results are derived from the comparison of an analytic criterion for gap-opening 
against SPH numerical simulations. In the simulations presented in section §3, the disk 
is represented by a collection of 2 x 10 5 SPH particles. In order to check that we have the 
resolution required and that our results do not depend on that, we run four additional 
simulations with different numbers of SPH particles; 5 x 10 5 , 8 x 10 5 and two with 10 6 . 

Although, the evolution of the orbital parameters of the binary are not exactly the 
same for high and low resolution simulations, the global structure of the disk remain 
practically almost unchanged. The region of parameters that the low and high resolution 
simulations covers, in the velocity spaces of parameters, is very similar (figure [TJ) and 
the classification of these simulations (as closed or opened) remain the same. For the 
initial condition that we run with 2 x 10 2 , 5 x 10 5 , 8 x 10 5 and 10 6 we find that the 
values of the velocity parameters converge with resolution (figure 12} • 




Figure 1. The cubic ratio between the sound speed of the gas and the rotational 
velocity of the binary-disk system, (c s /v) 2 plotted against the quadratic ratio between 
the rotational velocity of the isolated binary and the rotational velocity of the binary- 
disk system (ubm/f) 2 - The black dots are runs with 2 x 10 5 SPH particles and the 
black open circles are runs with 1 x 10 6 SPH particles. Also we plot all our others 
simulations with red {opened simulations) and blue {closed simulations) points for 
comparison. The green line is the interface between this two populations of parameters 
that is predicted by our analytic gap-opening criterion. The slope of the interface is 
the function /(A0, ao, (3q 1 a as ). The two simulations above (below) the green curve 
with 2 x 10 5 and 1 x 10 6 SPH particles are simulations with the same initial condition. 
Although, the evolution of the orbital parameters of the binary are not exactly the same 
for high and low resolution simulations, the global structure of the disk remain almost 
unchanged and the region of parameters that the low and high resolution simulations 
covers is very similar. The classification of these simulations (as closed or opened) 
remain the same for high and low resolution. 
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Figure 2. The cubic ratio between the sound speed of the gas and the rotational 
velocity of the binary-disk system, (c s /v) 2 plotted against the quadratic ratio between 
the rotational velocity of the isolated binary and the rotational velocity of the binary- 
disk system (v^-m/v) 2 . In this figure we plot four simulations with the same initial 
condition but different number of SPH particles: 2 x 10 5 , 5 x 10 5 , 8 x 10 5 and 1 x 10 6 . 
The green line is the interface between the simulations where a gap is opened in the 
disk (opened) and the simulations where the disk dont have a gap (closed). Although, 
the evolution of the orbital parameters of the binary are not exactly the same for 
high and low resolution simulations, the global structure of the disk remain almost 
unchanged and the region of parameters that the low and high resolution simulations 
covers is very similar. The classification of these simulations opened remain the same 
for high and low resolution. 



For the high resolution simulations we find a small shift to higher values of (c s /v) 3 . 
This shift can change the value of the slope of the threshold between the closed and the 
opened simulations (/) in a %10 or less. Therefore, the overall conclusions for the testing 
of our analytic gap-opening criterion against SPH simulations remain unchanged. 
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